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The change of the Josephson supercurrent density j s of 
a weakly-connected granular superconductor in response to 
externally applied arbitrary thermal gradient VT (nonlinear 
Seebeck effect) is considered within a model of 3D Joseph- 
son junction arrays. For VT > (VT) C , where (VT% is es- 
timated to be of the order of ~ 10 4 A"/m for YBCO ceram- 
ics with an average grain's size d ~ 10/wn, the weak-links- 
dominated thermopower S is predicted to become strongly 
VT-dependent . 



A linear Seebeck effect, observed in conventional and 
high-T c ceramic superconductors (HTS) n and attributed 
to their weak-links structure (see, e.g.,Er0 and further 
references therein), is based on the well-known fact that 
in a Josephson junction (JJ) the superconducting phase 
difference A(f> depends only on the supercurrent density 
j s (according to the Josephson relation j s — j c sin Afa 
where j c is the critical current density). When a small 
enough temperature gradient VT is applied to such a 
J J (with the normal resistivity p n ), the entropy-carrying 
normal current with density j n = SoVT/ p n is generated 
through such a junction, where So is the thermopower (a 
linear Seebeck coefficient). This normal current density 
is locally cancelled by a counterflow of supercurrent with 
density j s , so that the total current density through the 
junction j = j n + j s =0. As a resultu, the supercurrent 
density j s — —j n generates a nonzero phase difference 
Acj) via a transient Seebeck thermoelectric field Et — 
Pnjn = S'oVT induced by the temperature gradient VT. 
If in addition, an external current of density j e also passes 
through the weak link, a non-zero voltage will appear 
when the total current density exceeds j c , i.e., for j e — 
j c ± S VT/p„. 

In the present Letter, using a zero-temperature 3D 
model of Josephson junction arrays, a nonlinear analog of 
the thermoelectric effect (characterized by a non-trivial 
VT-dependence of the Seebeck coefficient S) in granular 
superconductors is considered. The experimental condi- 
tions under which the predicted behavior of thermopower 
can be observed in YBCO ceramics are discussed. 

The so-called 3D model of Josephson junction arrays 
(which is often used to simulate a thermodynamic be- 



havior of a real granular superconductor) is based on the 
well-known tunneling Hamiltonian (see, e.g.,uE3) 
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[1 - COS <f>ij(t)], 



(1) 



and describes a short-range interaction between N su- 
perconducting grains (with the gauge invariant phase 
difference faj(t), see below), arranged in a 3D lattice 
with coordinates n = (xi,yi,Zi). The grains are sep- 
arated by insulating boundaries producing Josephson 
coupling Jij which is assumedcl to vary exponentially 
with the distance r\j between neighboring grains, i.e., 
Jij(rij) = J(T)e~ K ' Tij . For periodic and isotropic ar- 
rangement of identical grains (with spacing d between 
the centers of adjacent grains), we have k = (4, 4, 4). 
Thus d is of the order of an average grain (or junction) 
size. 

In general, the gauge invariant phase difference is de- 
fined as follows 

4>ij{t) = <kj(0) - Aij(t), (2) 

where <pij(0) = (pi — 4>j with fa being the phase of the 
superconducting order parameter, and Aij(t) is the so- 
called frustration parameter, defined as 



A i3 {t) 



2ir 
$ 



A(r, t) ■ dl, 



(3) 



with A(f, t) the (space-time dependent) electromagnetic 
vector potential; $o — h/2e is the quantum of flux, with 
h Planck's constant, and e the electronic charge. 

As is knownE3il3, a constant electric field E applied to a 
single JJ causes a time evolution of the phase difference. 
In this particular case Eq.(2) reads faj(t) = faj(Q) + 
u/ij(E)t where Uij(E) — 2eE ■ r^/fi with f\j = fi — fj 
being the distance between grains. If, in addition to the 
external electric field E, the network of superconducting 
grains is under the influence of an applied magnetic field 
H, the frustration parameter Aij(t) in Eq.(3) takes the 
following form 
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Here, Rij = (r*i + fj)/2, and we have used the conven- 
tional relationship between the vector potential A and 
(i) constant magnetic field H = rot A (with dH/dt = 0) 
and (ii) homogeneous electric field E = —dA/dt (with 
rotE = 0). 

There are at least two ways to incorporate a thermal 
gradient VT dependence into the above model. Namely, 
we can either invoke an analogy with the above-discussed 
influence of an applied electric field on the system of 
weakly-coupled superconducting grains or assume a di- 
rect VT dependence of the phase difference (as it was 
recently suggested by Guttman et altil). For simplicity, 
in what follows we choose the first possibility and as- 
sume an analogy with the conventional Seebeck effect. 
Recall that application of a temperature gradient VT 
to a granular sample is known to produce a thermoelec- 
tric fieldEm Et — SqVT, where So is the so-called linear 
(VT- independent) Seebeck coefficient. Assuming that in 
Eq.(4) E = Et, we arrive at the following change of the 
junction phase difference under the influence of an ap- 
plied thermal gradient VT 
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with the frustration parameter 



A 3 {t) = —{H AR tj ) 
<p 



265*0 ^ 
t-^T ■ njt. (6) 



As we see, the above equation explicitly introduces a di- 
rect VT dependence into the phase difference, expressing 
thus the main feature of thft-sp-called thermophase effect 
suggested by Guttman et alia. Physically, it means that 
the macroscopic normal thermoelectric voltage V cou- 
ples to the phase difference on the junction through the 
quantum-mechanical Josephson relation V cx dA<j)/dt. 
Later on we shall obtain a rather simple connection be- 
tween the thermophase coefficient St = dAcf>/dAT and 
the conventional linear Seebeck coefficient So- 

To consider a nonlinear analog of the Seebeck effect 
(characterized by a VT-dependent thermopower S), we 
recallll3'EJ that within the model under consideration the 
supercurrent density operator j s is related to the pair po- 
larization operator p as follows (V is a sample's volume) 



f 1 dp 1 
where the polarization operator itself reads 



p = 



N 

E 



(7) 
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Here q\ = —2eni with ni the pair number operator, is 
the coordinate of the center of the grain. 

Finally, in view of Eqs.(l)-(8), and taking into account 
a usual " phase-number" commutation relation, [cf>i,nj] = 
iSij, we find 
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for the thermal gradient induced value of the averaged 
supercurrent density. Here a temporal averaging (with a 
characteristic time r) accounts for a change of the phase 
coherence during tunneling of Cooper pairs through the 
barrier, while integration over the relative grain positions 
fij is performed bearing in mind a short-range charac- 
ter of the Josephson coupling energy, viz. Jiji^ij) = 
J{T)f(x l3 )f{ VlJ )f{ ZlJ ) with f(u) = e-"/ d . 

To discuss a true VT induced thermophase effect only, 
in what follows we completely ignore the effects due to 
a nonzero applied magnetic field (by putting H = in 
Eq.(6)) as well as rather important-in granular supercon- 
ductors " self- field" effects (see OL3 for discussion of this 
problem) and assume that in equilibrium (initial) state 
(with VT = 0) < j s >= 0, implying thus </> y (0) = 0. 
The latter condition in fact coincides with a current,den- 
sity conservation requirement at zero temperatures. As 
a result, we find that an arbitrary temperature gradi- 
ent V K T = AT/ Ax, applied along the x-axis to the 
Josephson junction network, induces the appearance of 
the corresponding (nonlinear) longitudinal supercurrent 
with density 



j s (AT) ^< f B <y x T) >= j G(AT/AT ), (10) 



where 



G(z) = 



l + z 2 



(11) 



Here, j = 2eJNd/hV, (V a ,T) = AT /A.t = H/2edTS , 
and z = AT/AT . 

Notice that for a small enough temperature gradient 
(when Va;T <c V^Tj), we recover a conventional lin- 
ear Seebeck dependence j s (V x T) = a(T)SoV x T with 
a(T) = {2ed/h) 2 J(T)NT/V. On the other hand, for 
this result to be consistent with the above-discussed con- 
ventional expression j s = SqV x T '/ p n , zero temperature 
coefficient a(0) should be simply related to the specific 
resistance p n . Let us show that this is indeed the case. 
Using J(O) = hAo/4e 2 R n for zero-temperature Joseph- 
son energy (where Ao is a zero-temperature gap param- 
eter), V ~ Nld 2 for sample's volume (with I being a 
relevant sample's size), and taking into account that the 
normal state resistance between grains R n is related to p n 
as follows, p n ~ (d 2 /l)R n , the self-consistency condition 
a(0) = 1/ p n yields t ~ {I / d) 2 {Ti / A ) for the characteris- 
tic Josephson time. 

As it follows from Eq.(10), above some threshold value 
(V X T) C ~ 0.25(Va:T)o the supercurrent density starts 
to substantially deviate from a linear law suggesting 
thus the appearance of nonlinear Seebeck effect with 
VT-dependent coefficient S{V X T) = S /(l + z 2 ) where 
z = V X T/(V X T) and So = S(0). Let us estimate an 
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order of magnitude of this threshold value of the ther- 
mal gradient needed to observe the predicted nonlin- 
ear behavior of the thermopower in weak-links-bearing 
HTS. Using So ~ 0.5/j.V/K and A /k B ^ 90K for 
thermopower and zero-temperature gap parameter in 
YBCO, and I ~ 0.5mm for a typical sample's sizeau, we 
get r ~ UUjLs for the characteristic Josephson tunnel- 
ing time (CfB), and (V X T) C ~ 10 4 K/m for the thresh- 
old thermal gradient in a granular sample with an av- 
erage grain (or junction) size d ~ lOfim. Besides, tak- 
ing J(0) ~ A for a zero-temperature Josephson energy 
(in samples with R n ~ h/Ae 2 ), we arrive at the follow- 
ing reasonable estimate of the weak-links-dominated crit- 
ical current density j — 2eJ/hld ~ 10 3 A/m 2 in YBCO 
ceramics. We believe that the above estimates suggest 
quite an optimistic possibility to observe the discussed 
nonlinear behavior of the thermoelectric power in (real 
or artificially prepared) granular HTS materials and hope 
that the effects predicted in the present paper will be 
challenged by experimentalists. 

In conclusion, let us obtain the connection between 
the conventional (linear) Seebeck effect and the above- 
mentioned thermophase effect (which is. linear by defi- 
nition). According to Guttman et alliJ, the latter ef- 
fect is characterized by a nonzero transport coefficient 
St = dAtfi/dAT. In our particular case (with (f>ij(0) = 
and H = 0), it follows from Eqs.(5) and (6) that 




Hence, within our approach the above two VT induced 
linear effects (characterized by the transport coefficients 
St and So, respectively) are related to each other as fol- 
lows 




To summarize, the change of the Josephson supercur- 
rent density of a granular superconductor under the in- 
fluence of an arbitrary thermal gradient (a nonlinear See- 
beck effect) was considered within a model of 3D Joseph- 
son junction arrays. A possibility of experimental obser- 
vation of the predicted effect in HTS ceramics was dis- 
cussed. 
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